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2 ' Abstract 

The isospin-breaking correlator of the product of flavor octet vector currents, 



n 38 (g 2 ) = ifd 4 x exp(iq.x) < 0\T(V£(x)V*(0))\0 >, is computed to next-to- 
next-to-leading (two-loop) order in Chiral Perturbation Theory. Large cor- 



rections to both the magnitude and (^-dependence of the one-loop result are 
found, and the reasons for the slow convergence of the chiral series for the 
correlator given. The two-loop expression involves a single 0(q e ) countert- 
erm, present also in the two- loop expressions for II 33 (g 2 ) and II 88 (g 2 ), which 
counterterm contributes a constant to the scalar correlator II 38 ^ 2 ), defined 
by II 38 (q 2 ) = {q^qu — q 2 g^ u )IL 3S (q 2 ). The feasibility of extracting the value 
of this counterterm from other sources is discussed. Analysis of the slope of 
the correlator with respect to q 2 using QCD sum rules is shown to suggest 
that, even to two-loop order, the chiral series for the correlator may not yet 
be well-converged. 
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I. INTRODUCTION 



In the last decade, following the appearance of the classic papers of Gasser and Leutwyler 
i Ij-Ql , numerous treatments of low-energy hadronic properties employing the methods of Chi- 
ral Perturbation Theory (ChPT) have appeared (for an excellent recent review, see Ref. M). 
In the bulk of these treatments, the chiral expansion has been carried out to next-to-leading 
(one-loop) order (i.e. 0(q 4 ) in the usual chiral counting). Expressions for hadronic observ- 
ables, to this order, incorporate the constraints of current algebra and, in addition, provide 
the leading corrections to these constraints in a transparent and unambiguous manner. For 
many processes (see again Ref. {§) corrections to leading order results are ~ 20 — 30%, and 
truncating the full chiral series to this order, in consequence, appears well-justified. This is, 
however, not universally the case. For example, the one-loop amplitude for 77 — > 7r°7r° [|5]J|, 
which vanishes at leading order, differs significantly from experiment even near threshold. 
The same is true of the spectral function of the vector current correlator IT^, (q 2 ), where 

K (? 2 ) =ij*x exp(ig.x) < 0\T(V«(x)V u b (0))\0 >= (q,q u - q 2 9flu )U ab (q 2 ) (1.1) 

with the standard flavor octet vector current, = q^l^q, again even rather near nn 
threshold J/J . Even more dramatic is the case of the process i] — > 7r°77, for which the 
predicted one- loop branching ratio is a factor of ~ 170 smaller than the Particle Data 
Group [|l(| value. In the first two cases, the discrepancies between the one-loop results and 
experiment are a result of the fact that the leading order contributions vanish. Corrections 
at 0(q 6 ) are not unexpectedly large, and recent calculations to two-loop order, by Bellucci 
et al. [|TTJ for 77 — ► 7r°7r°, and by Golowich and Kambor for the spectral functions of 
11^ and II®®, demonstrate that inclusion of the 0(q 6 ) corrections to the 0(q 4 ) one-loop 
results brings the theoretical predictions nicely into accord with experiment for q 2 less than 
~ (8 — 9) m 2 . The importance of two- loop contributions, even rather near threshold, has also 



been demonstrated for the photon vacuum polarization function in Ref. ||12|| . The situation 
for r] — > 7r°77 (which most closely resembles the case at hand) will be discussed in more 
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detail below. Other examples of the necessity of including 0(q 6 ) contributions, in the odd 
intrinsic parity sector of C e ff, are also known, specifically 7r° — > 77*, r\ — > 77* pl3| , [14|j and 

77T + — > 7T + 7T , 7] — > / -f7T + 7T~ [0-|l5[ . 



In the present paper we will study the convergence of the isospin-breaking vector current 
correlator, H 3 ^ u , to two- loop order. We will show that, as for the chiral series of the ampli- 
tude for r] — > 7r°77, that of the correlator, IT^,, is quite poorly converged to one-loop order, 
and discuss the physical reasons for this similarity. We will also discuss evidence that, even 
to two-loop order, the latter series is not yet well-converged. It should be stressed that the 
correlator in question is of interest not only as an example of a quantity for which the chiral 
series is slowly converging, but is also of relevance to ongoing debates concerning the role 
of isospin-mixed vector meson exchange in isospin-breaking and, more particularly, charge- 
symmetry-breaking, observables in few-body systems (see Ref. [|16| for a discussion of a 
number of the contentious issues and list of other relevant references). Here the point is that 
one may choose the vector currents, rescaled by gv/^v (where my, gv are the correspond- 
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ing vector meson masses and decay constants, the latter defined via < 0|V r A f|V Aa (A) >= -^7^ 
with a the flavor and A the polarization label of the vector meson) as interpolating fields for 
the vector mesons. The isospin-breaking correlators Ilj^ and II™ then provide information 
on the g 2 -dependence of the off-diagonal elements of the vector meson propagator matrix, 
for this choice of interpolating fields. While the off-shell behavior of such propagator matrix 
elements is, in general, interpolating-field-dependent, one could couple the results for the 
propagator to those for the corresponding nucleon-vector meson vertices, obtained using the 
same choice of vector meson interpolating fields, to produce the relevant isospin-breaking 
contributions to NN scattering S-matrix elements, such S-matrix elements being indepen- 
dent of the choice of interpolating fields Finally, it should also be pointed out that 
the spectral function of 11^ is, at least in principle, measurable experimentally, though the 
accuracy required to extract it makes this a rather moot point, at present. The possibility 
of this extraction rests on the observation that the isovector vector current matrix elements 
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< Tr^V^O > (k = 1,2,3) receive isospin-breaking contributions only at second order in 
{ma — m u ) ]T^,[n|, whereas < 7T7t|V^|0 > is non-zero already at Oim^ — m u ). This means 
that the deviation of the ratio of vector spectral functions measured in t~ — > v t tt~7t and 
e + e~ — > 7T + 7r _ from that predicted by isospin symmetry is (up to corrections for the heavy 
quark pieces of the electromagnetic (EM) current) a direct measure of the spectral function 
of 11^. Since these effects will be seen to be of order a few times 10 -4 they are, however, 
well outside the reach of current experiments, for which cross-sections below the resonace 
region in e + e _ — > tt + tt~ are known typically to an accuracy of only ~ 10%. 

The remainder of the paper is organized as follows. In Section II we record the relevant 
terms of the effective Lagrangian at orders 2, 4 and 6 in the chiral expansion, and discuss 
the general, diagrammatic structure of the one-loop and two-loop results. In Section III we 
describe briefly some details of the calculations and quote the full one- and two-loop results 
for the contributions identified in Section II. Detailed formulae for the loop integrals entering 
these expressions are relegated to the Appendix. Since these integrals have been discussed 
in considerable detail elsewhere (see, for example Refs. |7],[2tJ), the Appendix will be rather 
brief, and the reader is referred to the references just cited for further details. Section IV 
provides a discussion of the results, in particular the physical origin of the slow convergence 
of the chiral series for the correlator to one-loop. In Section V, a rough estimate of the single 
0(q 6 ) counterterm appearing in the 0(q 6 ) corrections to the one-loop result is given, based 
on a QCD sum rule analysis of the correlator, and the possibility of independent estimates 
of this low-energy constant from other sources discussed. The issue of the convergence of 
the chiral series to two-loop order is also treated. Finally, in Section VI, we summarize our 
conclusions. 
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II. THE CHIRAL LAGRANGIAN TO 0(Q 6 ) AND THE STRUCTURE OF THE 
CONTRIBUTIONS TO THE CORRELATOR 



The leading terms in the low-energy, chiral expansion of the correlator, IT^,, may be 
obtained from the effective chiral Lagrangian, £ e //, which may be written in the form 

C eff = C^+C^ + C^ + --- (2.1) 

where the superscripts denote the chiral order. The general form of and in the 
presence of external scalar, pseudoscalar, vector and axial vector sources, is given in Ref. 
|TJ . Since we are interested only in the correlator of vector currents we may set the external 
pseudoscalar and axial sources to zero and the external scalar source to 2Bq M, where M is 
the current quark mass matrix and Bq the usual parameter, appearing in and related 
to the value of the quark condensate. One then has, explicitly, for and Q 

£ (2) = fTriD^UD^tf) + -f 2 Ti[B M(U + U j )] (2.2) 
4 2 

and 

£ (4) = L^TriD^UD^)} 2 

+ L 2 Tr(D l ,UD^) r Tr(D f "UD u U^ + L 3 Tr(D^U ] D^U D U U ] ' D V U) 

+ L4Tr(D D^U^ )Tr[2B M (U + U ] )} + L 5 Tr[2B (MU + U ] M)D^U ] D^U] 

+ L 6 [Tr[2B M{U + U^)}} 2 + L 7 [Ti[2B M(U - U^)]f 

+ L 8 Tr[AB 2 (MUMU + MU ] ' MU ] )} - iL 9 Tr[F fiv D"UD v U^ + F^ U D^U ] D U U] 
+ L l0 Tr[tfF^UF^] + H^F^F^ + F^F^] + H 2 Tr[AB 2 M 2 } . (2.3) 

In Eqns. (2.2) and (2.3), Bo is a mass scale related to the value of the quark condensate 
in the chiral limit, U = exp(zA • Tt/f) (with A the usual £77(3) Gell-Mann matrices and 7? 
the octet of pseudoscalar (pseudo-) Goldstone boson fields), / is a dimensionful constant, 
equal to f n in leading order, M is the current quark mass matrix, and _D M is the covariant 
derivative which, in the absence of external axial vector sources, takes the form 

5 



D^U = d^U -i[v^U}. (2.4) 

The vector field strength tensor, F^ u , occuring in Eqn. (2.3) is defined by F^ u = d^v v — 
dvV^ — i[Vfj,,v u ], where t> M = ju* with the octet of external 577(3) vector fields. For 
the case at hand we require only the external sources and vf, and hence the last term 
in F^ u vanishes. Note that one would have to supplement Eqn. (2.3) with additional terms 
involving Tr^F^) if one wished to treat the correlator IT^, but, as these terms do not enter 
the calculation of IT^,, we have not explicitly displayed them in (2.4). Note also that, in 
writing the form (2.3) for additional terms which vanish as a consequence of the lowest 
order equation of motion have been omitted. In performing calculations to 0(q 6 ) one would, 
in general, also have to include these terms. However, since the effect of their presence can 
always be absorbed into a redefinition of the coefficients occuring in , we may drop 

these terms from the outset. 

The general form of in the presence of external sources, has been determined recently 



by Fearing and Scherer f22]| . The full expression, however, contains 111 terms of even intrinsic 
parity and 32 of odd intrinsic parity, and so will not be recorded here. In fact, only one 
combination of the terms from actually enters the present calculation. It is easy to see 
why this is the case. As pointed out in Ref. J7|, since to 0(q 6 ) only terms from zeroth 
order in the meson fields contribute to vacuum expectations, and since only four terms of the 
143 mentioned above contain terms zeroth order in the meson fields and second order in the 
external vector sources and vf,, only these four terms can contribute to the correlators 11^, 
with a, b = 3, 8. thus reduces, for the purpose of computing such vacuum correlators to 
0(q 6 ), to § 



K x Tr{D x f^D x U, v ) + K 2 Tr{D»f>?D x f +Xu ) 
f K 3 Tr{f +flu f!?x+) + K 4 Tr(/ + ^/f)Tr( X+ ) (2.5) 



where 



X+ = 2B u{U ] M + MU)u ] . (2.6) 

In Eqns. (2.5), (2.6), u = exp(i\ ■ 7r/2/) is the usual square root of the matrix U defined 
above, the covariant derivative of /+" reduces to 

Dif? = d x fF-i[v x ,f!?\ (2.7) 

in the absence of external axial vector sources, and all other notation is as defined before. 
To zeroth order in the meson fields, f+* and x+ are equal to 2v fiu and ABqM, respectively. 
The terms involving K\, K% and K± obviously contain no pieces involving both and vf, 
and hence do not contribute to the correlator If^,. Only the K s term survives. To facilitate 
comparison with Ref. |7j] we introduce the rescaled version of the low-energy constant (LEC) 
K 3 , Q = 4K 3 . 

One may now easily characterize the full set of contributions to the correlator nj^, to 
0(q 6 ). Generically these are of two types, corresponding to the two ways in which terms 
involving the product vfp® can arise in the expansion of exp(i / cftxCefflv^]): (1) those terms 
arising from the second order term in the expansion of the exponential and hence generated 
by pieces of £ e // first order in the external vector sources, and (2) contact terms, arising 
from the first order term in the expansion of the exponential, and hence generated by those 
pieces of C e ff second order in the sources. The resulting contributions to Ilj^ are depicted 
graphically in Figures 1-6. In these figures the left-hand current line carries momentum, 
flavor and Lorentz indices q, 3 and fi, and the right-hand current line, similarly, the indices 
q, 8 and v. Open circles enclosing a cross appearing in Figs. 1-5 denote those vertices 
generated by and the open box enclosing a cross in Fig. 6 the vertex (proportional to 
Q) generated by C^ 6 \ All other vertices are understood to be from C^ 2 \ Fig. 1 contains the 
full set of contributions of 0(q A ), Figs. 2-6 those of 0(q 6 ). Figs. 2 and 3 can be interpreted 
as dressing the internal propagators of Figs. 1(a) and 1(b). Additional graphs of the form 
4(a) and 4(b), in which the structures at the left- and right-hand current vertices have 
been interchanged have not been shown explicitly, but are understood to be present. Since 
the various vertices appearing in the figures can be read off from the expressions for C^ 2 \ 
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£( 4 ) and £( 6 \ it is a straightforward exercise in Feynmann diagrammatics to evaluate the 
correlator. Results for the various contributions depicted in the figures are presented in the 
next section. 

III. THE CORRELATOR IL™ TO ONE- AND TWO-LOOP ORDER 

In this section we record the results for the various contributions to the correlator 11^,, 
together with a few salient features of the calculations. All loop integrals required have 
been performed using dimensional regularization and can be expressed in terms of the basic 
integrals 

A(m 2 ) = [ 1 9 (3.1) 

v ; J (2-n) d k 2 -m 2 y ' 

and 

B[m\ q 2 ) = jf 1 log(l - q 2 x(l - x)/m 2 ) , (3.2) 

which are given explicitly in the Appendix. The auxiliary tensors T^ u (m 2 ,q 2 ) and 
Djj, v {m 2 , q 2 ), which occur frequently in the calculations are also described there. In what 
follows, the tensor decomposition 

T^{m 2 , q 2 ) = T.im 2 , q 2 ){q^q u - q 2 g^) + 2A(m 2 )g^ (3.3) 

and the relation 

^(m 2 ,g 2 ) = ^(m 2 ,g 2 ), (3.4) 

which follows from the expressions given in the Appendix, have been used to reduce the 
results to compact forms involving the integrals A, B and T±. The expression for T\ in 
terms of A and B may also be found in the Appendix. In order to streamline notation we 
will write T\(P) for T^mj^q 2 ), B(P) for B(m 2 P ,q 2 ) and A{P) for A{m 2 P ) in what follows, 
where P = K + , K~ , ir and rj and the masses are understood to be those given by the 



8 



leading order chiral relations m 2 = 2B m, m 2 K+ = B Q (m s + m u ), m 2 K0 = B (m s + ma) and 
m 2 = 2B (2m s + rh)/3, where m = (m u + m^)/2. 

Let us begin with the one-loop, C^g 4 ), result generated by the diagrams of Fig. 1. One 
may easily verify that there are no contributions to IT^, of the type 1(c), and that the 0(q i ) 
contribution is generated solely by differences between K + and K° loops of types 1(a) and 
1(b). The contributions of type 1(a) are 



T^(K + )-T^(K°) 



(3.5) 



and those of type 1(b) 



A(K°) - A(K A 



9^u ■ 



(3.6) 



The sum of these contributions yields the full 0(q A ) result for the correlator, 



W) 



(4) i\f2>, 



(3.7) 



which has, of course, the transverse structure required of the vacuum value of the covariant 
time ordered product of the conserved vector currents and V^, and can also be seen to 
be finite and manifestly independent of the renormalization scale, fi, from the form of T\ 
given in the Appendix. 

Turning to the contributions of 0(q 6 ), we begin with the insertion graphs of Figs. 2,3. 
These are known in terms of the one-loop contributions to the wavefunction renormalization 
constants and mass shifts of the internal (K + and K°) lines. The resulting expressions are 
considerably simplified if we include also the contributions of type 4(a) and 5(b) involving 
the LEC's L 4 and L 5 , since those contributions exactly cancel the terms involving explicit 
factors of L 4 and L 5 arising from Figs. 2(a) and 3(a). The resulting contributions to LT^ 
are then 



{q^qu-q 9iiv) 



2AP 



3 (t^K ) - T^K + )) (A(K°) + A(K + ) + A(ir) + A( v )) 



+ (T^K ) + T^K*)) (A(K°) - A{K + ) + 2^0 (A( V ) - A(ir 



iV3 



(5m 2 K+ )i^i oop 



($m 2 K o)i-i oop 



B(K C 



(A{K°) - A(K + )) (AA(K + ) + AA{K°) + 3A{n) + 3A( V ) 



+2V39 (A(K°) + A{K + )) {A{rj) - A(n)) 



(3.8) 



where 9 = \/3(md — m u )/A(m s — m) is the leading order tc-i] mixing angle and (5m 2 K+ )i-i oop , 
i— loop &re the one-loop corrections to the leading order K~^~ and squared masses, 
the expressions for which may be found in Ref. [jj]. The remaining terms of types 4(a) and 
5(b), which involve the LEC's L 910 , then yield the contribution 



{%<lv - q g^u) 



24/ 2 



2AiL 9 q 2 (T X (K°) - T 1 (K + )) - 48* (L 9 + L 10 ) (A(K°) - A(K^ 



(3.9) 



The remaining contributions to Ilj^ are: (1), from 4(b) 

r./5i r . 



\V3' 




[24P_ 


[( 
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T^K ) - T l (K + )j (9A(K°) + 9A(K+) + 9A(ir) + 3A(r)) 
+ (T\(K°) + T 1 {K + )) (7A{K°) - 7A{K + ) + 2V39 (A( V ) - A{n)) 
+6Ti(tt) (A(K°) - A(K + 



fJbV 



24/ 2 



(A(K°) - A(K + )) (32A(K Q ) + 32A(K + ) + 3(L4(vr) + 6A(rj)) 



+AV30O (A(K°) + A(K + )) (A( V ) - A(tt)) 



(3.10) 

(2), from 4(c), (where, owing to the structure of the loop integrals, only contributions with 
the central vertex from the kinetic portion of survive) 

{%q v - q 2 9nv) 



r V5- 




[24/ 2 _ 





3q 2 (Ti(^°) - T 1 {K+)) (^(tt) + Tx(^ ) + T^K 
-6 (^(K ) - T 1 (K + )) (A(K°) + A(K + ) + A(n 
-6 (AiK ) - A(K + )) (^(tt) + T^K ) + T 1 {K+) 

12 (A(K°) - A(K + )) (A(tt) + A(K°) + A(K 



24/ 2 



(3.11) 
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(3), from 5(b), 



"24/2 



and (4), from Fig. 6, 



(A(K°) - A{K + )) (l6A(K°) + 16A(K + ) + 15A(n) + 3^(77)) 



+2^0 (A(K°) + A(K + )) (A(ri) - A(tt)) 



(3.12) 



(3.13) 



Adding the results of Eqns. (3.8) through (3.13), we obtain, for the full 0(q e ) contribution 



to Uf v , 



n^(g 2 )l (6) = (<v?,-<?V) 



iVS rfx 2 v B(K + ) , ..... 

{{(> m K+)l-loop — 2 {Om K o)l-loop 2 J 



+ 



AQ 

(m 2 K o - m 2 K +) - 48i(L 9 + L 10 ) 



B(K°) 

J 2 ' 



+q 



v/3/ 2 



24/2 



(A{K°) - A(K+)) 



24/ 2 



{Ti(K°) — Ti(K + )} 



x{3 (7\(7r) + T^K ) + T 1 (K + )) + 24tL 9 } 



(3.14) 



As expected, the non-transverse contributions appearing in (3.8), (3.10), (3.11) and (3.12) 
have all cancelled. In Eqn. (3.14), we may replace the lowest order expressions for the 
meson masses appearing in T x , B and A with the physical masses, to the order we are 
working. This is not, however, true of Eqn. (3.7). If we wish to combine the results of (3.7) 
with those of (3.13), we must re-express the leading-order squared masses occuring on the 
RHS of (3.7) as the differences of the corresponding one-loop expressions (which can then 
be set to the physical masses when working to 0(q 6 ) overall) and the 0(q A ) corrections, 
{5rn 2 K+ )-i-i oop , {5rn 2 K o)i-ioop- To 0(q 6 ) overall it is then sufficient to expand Ti (m 2 K+}0 , q 2 ) 
about the physical values of the squared masses, M^-+, , to first order in {5m 2 K+ )i-i oop , 
{5m 2 K o)i-i oop . The derivative of Ti(m 2 ,g 2 ) with respect to m 2 which is required here can 
be obtained from the expressions in the Appendix. The terms first order in {5m 2 K+ )i-i OCfp , 
(Sm 2 K o)i-ioo P which result turn out to cancel those in (3.13). 
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Before recording the final 0(q 4 )+0(q 6 ) result for H 3 ^ u , we must discuss the renormaliza- 
tion prescription implicit in Eqn. (3.14) (as pointed out above, (3.7) is already finite and 
scale- independent). The loop integrals T\(P), A(P) and the LEC's £9,10 and Q all contain 
divergences as d — > 4. Those of Z/9,10 are already known from the renormalization of C e ff to 
0(q A ) |TJ, and those of Ti(P), A(P) are given in the Appendix. Note that, since the vertices 
arising from and involving £9,10 appear in divergent loop graphs, one must go beyond 
the expressions for Lg^o used in 0(q 4 ) calculations and include the next terms in the Laurent 
expansions of these LEC's in terms of the variable A = -737 — \ (log(47r) — 7^ + 1) 



Lk = /- 4 [lg>Qi)\ + 4°V) + ^i" 1 ^- 1 ] , (3.15) 

where \x is the dimensional regularization renormalization scale and 7^ is Euler's constant, 
and, in the more familiar notation of Ref . |1| , 

4 1} = r fc and Lf=U k . (3.16) 

One would similarly require, in general, an expression for the 0(q 6 ) LEC, Q, of the form 



<f-4 



q = (y] [Q^^)y + q w (^ + Q (0) (ri\ , (3.17) 

in order to absorb all divergences in two loop calculations. From (3.14) it follows that 

Q( 2 V) = and Q w (fj) = 3 (4°V) + jj$>(ji)) , (3.18) 



in agreement with the results of Ref. 0, whose notation we have followed in the expressions 
above. Note that the LEC's L910 occur ! as claimed earlier, only at 0(q 6 ), and in the fixed 

combination, Q (o) (/-0 = Q (0) G") ~ 3 (4°V) + 4o (A*)), with the LEC Q {0) (v)- The 

scale dependence of the various LEC's is discussed in detail in Ref. and will not be 
repeated here. 

Given the expressions for the divergent pieces of the loop integrals Ti(P) and A(P) 
and the LEC's £9,10, one may now easily verify that the quantities enclosed in braces in 
Eqn. (3.14) contain no divergences as d — > 4. As such the 0[d — 4) terms in the expansions 
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of the Ti(P), B(P) and A(P) are not required, and hence have not been recorded explicitly 
in the Appendix. The result is finite after the renormalization of Q given in (3.17), (3.18) 
above. The results of Ref. for the scale-dependence of the LEC's also allows one to check 
that the 0(q 6 ) result, (3.14), is scale-independent. 

Using the expressions (3.16), (3.17) and (3.18) (with T 9 = —T w = 1/4 from Ref. jlj), the 
explicit form for T±(P) given in the Appendix, and rewriting (3.7) in terms of the physical 
K + and K masses as described above, one obtains the following compact form for the scalar 
correlator, II 38 (g 2 ), valid to 6 th order in the chiral expansion: 

-2iB{M 2 Kl q 2 



tf\q 2 ) = ^-{Ml -Ml + ) QCD 



q 2 



1 + 



2hf - i (B 21 (M 2 , q 2 ) + 2B 21 (M 2 K , q 2 )) - ^ log^M* ///) 



/r(0) , r(0K 1fi 



(3.19) 



where M\ is the average of the non-EM portion of the physical K + and K° squared masses, 
(Mj^o — Mk+)qcd is the non-EM contribution to the kaon mass-squared splitting, and the 
auxiliary quantity B 2 \(M 2 , q 2 ) is defined, in terms of B(M 2 ,q 2 ), by 

ftl (MV) = l(l-^)s ( M^V^. (3.20) 

For later reference we record here the expression for the imaginary part of Il^g 2 ), valid 
for q 2 < 4M|-, and to sixth order in the chiral expansion: 

Imn (g ) = 192^2/2 Re(^B(M^,g)J 

/ 4M 2 \ 3/2 

x^l-^j % 2 "4M 2 ). (3.21) 
This expression follows straightforwardly from (3.19) and the results quoted in the Appendix. 



IV. DISCUSSION OF THE RESULTS 

The expression (3.19) gives the full result for II 38 , (g 2 ), valid to sixth order in the chiral 
expansion. The functions iB(M 2 ,q 2 ), and hence also iB 2 i(M 2 ,q 2 ), have cuts beginning at 
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q 2 = AM 2 . Since q 2 = AM\ is, presumably, outside the range of validity of the chiral expan- 
sion, the imaginary part of II 38 is generated solely by the B 2 i(M 2 ,q 2 ) term in Eqn. (3.19), 
in the region of q 2 values of interest to us here. This contribution arises from graphs of the 
form 4(c) having 7r + 7r _ fields at the vertex and K + K~ or K°K° fields at the vertex. 
Because of the absence of a lowest order (C^) coupling of to n + ix~ , the one-loop result 
for ImIT 38 (g 2 ) is zero for q 2 < 4M|-; the tttt cut first enters only at two-loop order. 

As explained below, all quantities appearing in (3.19) are previously known, with the 
exception of the tree-level 0(q 6 ) LEC . Because of significant cancellation between the 
g 2 -independent part of the one- loop contribution and the term involving (L^ + L^), the 
q(o) term will undoubtedly play a significant role, particularly near q 2 = 0, and may even 
be large compared to these other constant terms individually. We will comment further 
on this question below, but for now will leave as unknown and investigate the size 
of the genuine 0(q 6 ) loop corrections to the one-loop result. In Figure 7 the real and 
imaginary parts of IT 38 (g 2 ) are displayed (less the constant contribution to ReIT 38 (g 2 )) 
for — 12 M 2 < q 2 < 12 M 2 , together with the one-loop result (which is purely real in this 
range). One should note that, although the results are displayed out to q 2 = 12 M 2 , the 
corresponding two- loop expression for II 33 (q 2 ) begins to deviate from experiment above 
q 2 ~ (8 — 9)M 2 fTj. Moreover, as will be discussed below, one should bear in mind that 
the range of validity of the two-loop expression for ImIT 38 (g 2 ) may not extend as far above 
threshold as does the that of the corresponding expression for ImII 33 (g 2 ). 

In arriving at the numerical results shown in Fig. 7, we have used the following input 
information. First, we follow standard practice in taking / ~ f n = 92.4 MeV. Second, 
having demonstrated the scale-independence of (3.19), we set \x = Mk, which simplifies the 
logarithmic terms, and use the values 4 0) (/i = M K ) = 0.0073 ± 0.0003 and = m k) = 

—0.0058 ± 0.0003 (see Refs. [0 and p^| ) (compatible with the resonance saturation result 
L^q , (// = m p ) = — |Lg ()Lt = m p ) [p4] , p5[ ). The errors on Lg°] are more significant in the 
combination + L±q, both due to the large cancellation between the central values and 
to the cancellation between the one-loop and genuine two-loop contributions, resulting in a 
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considerable variation in the precise value of II 38 (0). Indeed, including the uncertainties in 
Lg "* + Liq by allowing either LEC to vary within the quoted error bars, one finds the full 
one-plus-two-loop curve for the real part of the correlator in Fig. 7 is shifted up or down 
by ~ 0.5 x 10" 5 , with little change in shape (the imaginary part is not affected). Since, 
however, the contribution is expected to be dominant (see also the discussion below), 
this uncertainty is unlikely to be significant for the correlator as a whole and, as a result, 
we have used the central values to obtain the results shown in the figure. Finally, we have 
obtained the non-EM contribution to the kaon splitting, associated with (rrid — m u ) ^ 0, by 
subtracting the EM contribution, where the latter is evaluated as follows. In the past the 



EM subtraction has been made using Dashen's theorem [26 



(Mf + - Ml„) EM = (M| + - Ml) EM * (M| + - Ml) apt , (4.1) 

a result valid strictly only in the chiral limit. Recently, arguments have been advanced 
|27|-p9"| suggesting that the theorem receives significant corrections beyond leading order 
and we have, therefore, used the value 

(M^-M^) EM ^l.9(A^-M^ evt (4.2) 

suggested by the analyses of Refs. [|28 , 29| in arriving at {M 2 K+ - M 2 K0 ) QCD . 

One significant feature of the results is immediately obvious from Fig. 7: despite being 
higher order in the chiral expansion, the genuine loop contributions of 0(q 6 ) are actually, for 
most of the q 2 range displayed, even larger in magnitude than the one-loop, 0(q A ), result. 
Moreover, unlike the 0(q 4 ) result, which has a rather small variation with q 2 , the 0(q 6 ) 
corrections are strongly g 2 -dependent. We thus see that, independent of the value of , 
the chiral series for II 38 (q 2 ) is poorly converged to 0(q A ). 

The failure of the chiral series for II 38 (q 2 ) to be well-converged at one- loop order is 
actually not a surprise, given the similarity of the qualitative features of the 0(q 4 ) result, 
Eqn. (3.7), to those of the amplitude for rj — > 7r°77 MM- In the latter for II 38 (g 2 ), 

the leading, 0(q 2 ), contributions are zero, and the 0(q 4 ) LEC's, {L^}, do not contribute. 
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Moreover, loop contributions with internal n legs are suppressed by a factor of (ma — m u ), 
while loop contributions with internal K legs, which are not so suppressed, are instead 
suppressed by the natural smallness of the K loop integrals. The latter effect can be easily 
seen in the behavior of the loop integral function B(M 2 , q 2 ) near q 2 = 0: 

Thus, near q 2 = 0, the K loop integral is, e.g., suppressed by a factor of M 2 /Mx ~ 0.08 
relative to the corresponding ir loop integral. The result is that the one-loop prediction 
for the branching ratio of rj — > 7r°77 is a factor of ~ 170 smaller than that determined 
experimentally fL0 |. The reason for this discrepancy is well-understood, and is relevant to 



the case at hand. As is well-known |3(J,^4|,|25| , by making standard field choices, one may 
incorporate heavy resonances and the pseudoscalar octet into a single effective chiral La- 
grangian. Integrating out the heavy resonance fields then produces an effective Lagrangian 
for the pseudoscalars of the form given in Section II. The effect of the heavy resonances is to 
produce contributions to the LEC's. These contributions are fixed in terms of the parame- 
ters describing the couplings of the pseudoscalars and the heavy resonances in the original, 
extended, effective Lagrangian, these parameters, in turn, being fixed by comparison with 
experiment. One then finds that, where vector and axial vector resonances can contribute to 
a given their contributions practically saturate the observed values |24|,|25|]. Thus, for 
i] — > 7T°77, where the dominant contribution to the amplitude is known to be due to vector 
meson exchange [p|,|3T|, the absence of the in the one-loop result indicates the complete 



absence of the dominant contributions to the amplitude at this order (at least for the in- 
terpolating field choice for the vector mesons implicit in the standard construction). The 
effect of the vector mesons, in this case, first appears in the tree-level constants generated by 
terms from C^ 6 \ and these terms must, therefore, actually dominate the amplitude The 
situation for II 38 (g 2 ) is very similar. Here we again expect significant, probably dominant, 
vector meson exchange contributions, and the absence of the 0(q 4 ) LEC's from Eqn. (3.7) 
indicates that these contributions are not present in the 0(q A ) result. It is, therefore, likely 
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that the term in (3.19) will be the dominant one, at least at low q 2 , especially given the 
cancellation between the 0(q 4 ) result and the genuine loop corrections of 0(q 6 ) . (The 
term, of course, does not contribute to the slope of II 38 (q 2 ) with respect to q 2 , and this slope 
will, therefore, continue to be dominated by the the genuine two-loop 0(q 6 ) contributions.) 
We will discuss the possibilites for constraining in the next section. 

We close this section with a brief elaboration of our earlier comments on the experimental 
accessibility of the spectral function, p 38 (g 2 ), defined by 

lmU 38 (q 2 ) = 7r 6(q 2 - AM 2 ) p 38 (q 2 ) , (4.4) 

from a comparison of e + e~ — > 7r + 7r~ and r~ — ► z/ r 7r _ 7r data. The possibility rests on 
the presence of 0(m<f — m u ) isospin-breaking contributions in the former process, but their 
absence in the latter. As is well-known, the spectral function of the photon vacuum polariza- 
tion, p yi {q 2 ) 1 which is a linear combination of n 33 , II 38 , and II 88 , is directly proportional to 
the measured cross-section for e + e~ — > hadrons. Below q 2 = 9M 2 , only 7r + 7r~ intermediate 
states contribute. Since the coupling of V 8 to is already 0{md — m u ), one has, to 

0(m d - m u ), 

p"tf)=fPtf) + -j=fPtf), (45) 

and the deviation of p 11 {q 2 ) from the value expected based on t~ — > u T 7i^7T and isospin 
symmetry is, therefore, attributable completely to p 38 (q 2 )- The ratio r = 2p 38 (g 2 )/v / 3p 33 (g 2 ) 
thus represents the accuracy required in order to be able to extract p 38 {q 2 ) experimentally. 
The expression for p 38 (q 2 ) which follows from Eqns. (3.21) and (4.4) produces values of r 
of order 2 x 10~ 4 . It should be noted that the two-loop result for p 38 (q 2 ) is likely to be 
accurate only relatively close to threshold. This is because the graphs 4(c) which contribute 
to p 38 in this range of q 2 contain no rescattering. As is evident from the deviation of the 
one-loop result for p 33 from experiment (and hence from the full two-loop result) even at 
rather small q 2 (q 2 > 5M 2 ) - see Fig. 5 of Ref. [0 ) - such effects can be quite significant. If 
we use the p 33 results as a guide, such (yet higher order) effects might enhance p 38 (q 2 ) by a 
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factor of 2 or so in the vicinity of q 2 ~ 8 M 2 . This still leaves the isosp in-breaking correction 
factor, r, at ~ 4 x 10~ 4 , as mentioned in the Introduction, far outside the reach of current 
experiments. 

V. THE LEC <K°) AND THE CONVERGENCE OF THE CHIRAL SERIES TO 

TWO-LOOP ORDER 

The asymptotic behavior of the scalar correlator, II 38 (q 2 ) is known from the operator 



product expansion (OPE) to be, up to logarithmic corrections |32] 



U^fa 2 ) ~ ^ (m ^~ m » } (5 1) 

11 {q) -8n 2 {-q 2 ) ■ {bA) 

From (5.1), it follows that II 38 (g 2 ) satisfies the unsubtracted dispersion relation 

n 3 V) = r ds p38 2 (g) , . (5.2) 

J4M% s — q z — it 

As usual, this means that II 38 (q 2 ) and its derivatives with respect to q 2 at q 2 = can be 
written as negative moments of the spectral function p 38 (q 2 ), e.g., 

n 38 (0) = [°° ds (5.3) 

J4M2 s 



n 38 (0) = / ds . (5.4) 



An 38 (o) = r ds ^ s 

dq 2 UMl s 2 

The LHS's of these relations have chiral expansions which involve the quark masses and the 
LEC's appearing in C e ff. If one had experimental access to the spectral function, these 
relations (often called chiral sum rules) would serve to provide information on the LEC's. 
We have not written down the explicit form for these sum rules since, as pointed out above, 
p 38 {q 2 ) is unlikely to become experimentally available in the near future, but they are easily 
constructed from Eqn. (3.19). Eqn. (5.3), in particular, involves the new, unknown 0(q 6 ) 

LEC g (0) • 

Since we cannot realistically hope to constrain using (5.3), we must look for other 
ways to estimate its value. Perhaps the most favorable source of such an estimate would be 



the chiral sum rule, analogous to (5.3) above, for the difference of the II 33 and II 88 spectral 
functions, as derived in Ref. 0: 

' p 33 (s)-p 8S (s) \ 1 Ml 

\m-:\ s J 48tt 2 Ml 

*% log m max) + LfSi M K) ) + l%Mg,,,M K ) . (5.5) 



8tt 2 / 2 ° \Ml ) V y ^ iU v ~V 3/2 

This involves the 0(q 4 ) LEC's, and Li }, which are already rather well-known, and 
the at-present-unknown 0(q 6 ) LEC Q^. An analsyis of the sum rule (5.5) (in addition 
to other sum rules involving p 33 and p 88 ) is being performed by Golowich and Kambor, 
and should provide a useful estimate of , but at present this analysis has not been 
completed. One might be tempted to follow the path of estimating using the resonance 
saturation hypothesis, which has proven very successful for the 0(q A ) LEC's, and has also 
been employed in treating the 0(q 6 ) LEC's appearing in 77 — > 7r°7r° |H| and 77 — > 7r°77 
0. In the present case, however, the application of this method is more complicated than 
in previous situations since the term of interest in (involving the LEC Q) is generated 
only by graphs involving one 0(q 2 ) and one 0(q A ) vector meson vertex from the original, 
extended vector-plus-pseudoscalar effective Lagrangian. In order to fit the constants which 
determine the 0(q A ) vertices, one would have to do a detailed analysis of the vector meson 
EM decays which included the pseudoscalar loop corrections. While such an analysis would 
be of interest, given that the observed vector meson decay constant ratios show definite 
deviation from SU(3)p predictions, it is not available at present. We will, therefore, content 
ourselves with an alternate estimate based on a QCD sum rule analysis of the correlator in 
question. 

A sum rule analysis of the related correlator, II 3 ^, 

KM 2 ) = * J d 4 xex V (iq ■ x) < 0|T(V; 3 (x)V7(0))|0 > , (5.6) 
where Vjf = [u^ v u + d^ v d^ /6, has recently been performed |16| , updating the earlier anal- 



yses of Refs. |32| and |33| . Defining 
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n^(g 2 ) = i I d A xex P (iq ■ x) < 0|T(V; 3 (a;)^(0))|0 > , (5.7) 

with V$ ee s r y u s, and the scalar correlators U 3uJ (q 2 ), Il 3 ^(q 2 ) in analogy with II 38 (g 2 ), one 
then has 

n 38 (g 2 ) = y/3U^{q 2 ) - ^n 3 ^(g 2 ) . (5.8) 

a/3 



U 3lJJ (q 2 ) was analyzed in Ref. [16] by keeping terms of dimension six or less and working to 
first order in oem, &s and m q in the vacuum value of the OPE of the product of currents, 
V 3 (x)VJf(0). Including contributions to ImiT 3 ^ associated with the p, u, 0, p' and uj' 
mesons, one finds a very stable analysis for the correlator and one that, via the unsubtracted 
dispersion relation satisfied by II 3 ^, can be turned into a representation of the behavior of 
the correlator in the vicinity of q 2 = 0. Uncertainties in the values of the input four-quark 
condensates limit the numerical accuracy of this representation, but the values of n 3aJ (0) 
and ^n 3aJ (0) appear fixed, certainly to within a factor of 2 [1^]. If a similar analysis can 



be performed for n 3< ^(g 2 ), then we may use (5.8) to provide constraints on our two-loop 
representation of II 38 (q 2 ). 

The analysis of IT 3 ^(g 2 ) closely follows that of U 3w , so we will be rather brief here (the 
reader is referred Refs. [j32]-|3~4|| for technical details). It is immediately obvious that, owing 
to the flavor mismatch between the two currents, the dimension 2 and 4 contributions to the 
correlator are absent, at least up to and including terms of 0(a 2 ). The leading contributions 



are then of dimension 6 and, to 0{olem-, oc s ), have the following form |32| , |33 



< 



Q 6 



0| (u lalb \ a u - d lal5 \ a d) (s 7 Q 7 5 A a s) |0 > 

/ 2 1 - \ 

( y, < 0| [~-u-f a -f 5 u - -d-? a ^d) (s7 a 7 5 s) |0 > , (5.9) 

where Q 2 = —q 2 . The mixed flavor condensates appearing in (5.9) vanish in the standard 
vacuum saturation approximation and, being Zweig rule suppressed, are expected to be 
significantly smaller than analagous flavor diagonal four-quark condensates in any case. 
Estimates for such mixed flavor condensates were made in Ref. where, for example, it 
was found that, for q — d,s, 
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< 0| (gToTsA^j {<lTl^ a q) |0 > „ Q Q6 (5 1Q) 

< 0| (M7«75A a M) (M7 Q 7 5 A a M) |0 > 

As such, we should be able to safely neglect (5.9), which approximation then leads to 

IT 3 V) ~ \fm 3 "(q 2 ) . (5.11) 
The results of Ref. |16[ then imply 



n 38 (0) = (1.5 ± 0.4) x 10" 4 (5.12) 

-^n 38 (0) = (6.2 ± 2.2) x 10~ 4 GeV- 2 (5.13) 
dq 2 

where the quoted errors reflect uncertainties in the input values of the four-quark conden- 
sates. The overall scale of the results in (5.12) and (5.13) is set by the quark mass difference 
contribution to the kaon splitting, which has been determined from the experimental split- 
ting using the modified version of Dashen's theorem (Eqn. (4.2) above). Using the value 
determined, instead, by the unmodified form of Dashen's theorem would lower both numbers 
by ~ 20%. It should be stressed, for the sake of the discussion below, that attempting to 
raise the magnitude of the four-quark condensates sufficiently to lower the values in (5.12) 
and (5.13) beyond the lower bounds quoted there leads to instabilities in the analysis, i.e., 
the absence of a "stability window" for the extracted resonance parameters as a function of 
the Borel mass [16]; such values are, therefore, inconsistent, at least in the context of the 



analysis as presently performed. 

If we accept the approximations above, then (5.12) implies 

)(°) ~ (R fi -I- 9 71 y in" 5 



Q w ~ (8.6 ±2.7) x 10~ 5 , (5.14) 

where the errors relect both the uncertainties in (5.12) and those in the LEC's Lg '* and 
Ljq , though they are dominated by the former. As argued on physical grounds above, the 
q{o) term indeed almost completely dominates II 38 (0), being a factor of ~ 6 larger than 
both the 0(q 4 ) and 0(q 6 ) genuine loop contributions (which are comparable in magnitude, 
but opposite in sign). This, of course, raises questions about whether yet higher order 
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contributions are necessarily negligible. (Recall, e.g., that for i] — > 7r°77, the full vector 
meson dominance contribution to the partial width was a factor of 1.7 larger than that 
associated with only the 0(q 6 ) portions thereof ||.) If we use also the information contained 
in (5.13), then it, in fact, appears that there is good evidence for believing that higher order 
contributions must, indeed, be important. This statement follows from the observation 
that the slope of II 38 (g 2 ) with respect to q 2 is rather well-determined in (3.19), and is 
~ 9 x 10~ 5 GeV -2 , a factor of at least ~ 4 smaller than that given in (5.13). Small 
corrections to the dimension 6 contributions on the OPE side of the sum rule are incapable 
of altering this conclusion. The statement, to this order, is of course also independent of 
the value of Q^°\ which makes only a constant contribution to the scalar correlator. Thus, 
if (5.13) is even reasonably accurate, it clearly demonstrates that even to two-loop order 
the chiral series for the correlator is not yet well-converged. This is somewhat surprising, 
given the behavior of the flavor-diagonal correlators to two-loop order [0], but perhaps not as 
much so as one might, at first, think. Indeed, the slope of the correlator receives its largest 
contribution from the Lg term in (3.19), which contribution is associated with graphs of 
the type 4(a) in which the C^' vertex involves L 9 . Since the L 9 current vertices do not 
themselves break isospin, these graphs, like those which contribute to 0(q 4 ), involve only 
internal K lines and, as a result, the slope is suppressed by the smallness of the loop integral 
factor B(M 2 K ,q 2 ). As we saw for the correlator itself, in going from one- to two-loop order, 
when one creates the possibility of internal ir lines by going to higher chiral order, such 
"suppressed" contributions can be significantly enhanced. At present it is not clear whether 
or not this is actually the case here, but it is clearly worth further investigation. It will, 
in particular, be very interesting to compare the outcome of the analysis of the sum rule 
(5.5) with the estimate (5.14) for Q^°\ If the two agree, then one will be justified in having 
increased confidence in (5.13), as well as (5.12), and the case for the slow convergence of the 
chiral series for the correlator II 38 to two-loop order will be considerably strengthened. If 
not, it would point to some problem with the truncations usually made in applying the sum 
rule method, in the case of the correlator II 38 . 
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VI. SUMMARY 



In summary we have evaluated the mixed-isospin vector current correlator 11^ (g 2 ) to 
sixth order in the chiral expansion. The result is given in compact form in Eqn. (3.19), 
and involves the previously known 0(q 4 ) LEC's and , and a single combination of 
0{q % ) LEC's, Q( \ The results shows that (1) the genuine two-loop contributions to the 
correlator are, over much of the q 2 range considered, larger than the leading, 0(q 4 ) result 
and that (2) in contrast to the one-loop result, the two- loop expression has a very strong 
(^-dependence. An analysis of the correlator using QCD sum rules yields an estimate for the 
0(q % ) LEC, but at the same time indicates the likelihood that the two-loop expression 
for the correlator is not yet fully converged. Further work on the value of the LEC is 
required in order to clarify this issue. 
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APPENDIX: EXPLICIT EXPRESSIONS FOR THE VARIOUS LOOP 

INTEGRALS 



We list here the explicit forms of the various loop integrals which enter the results of 
Section III. A detailed discussion of most of the quantities listed below is given in Ref. EC 
and in Appendix A of Ref. @, to which the reader is referred for details. 
The scalar integral, A(m 2 ), already defined in the text, is given by 

m 2 log(m 2 /// 2 



A{m l 



d-4 



2m z X + 



16tt 2 



+ 0(d-4) 



(Al) 



with n the regularization scale and A as defined in the text. 

Defining the integrals B(m 2 , q 2 ), B^m 2 , q 2 ) and 5^(m 2 , q 2 ) associated with the graphs 
of Figs. 1(a) and 4(a-c) by 

cfik {1, k^, k^ky} 



{B, B^, B^ u } = J 



(2n) d [k 2 - m 2 ][{k - q) 2 - m 2 } ' 
one finds that the integrals B^ and B^„ occur only in the combinations 



(A.2) 



T M „(m 2 , q 2 ) = 45^(m 2 , q 2 ) - 2q^B v {m 2 , q 2 ) - 2q„B^m 2 , q 2 ) + q^q u B(m 2 , q 2 ) (A.3) 
D^m 2 , q 2 ) EE 2B^(m 2 , q 2 ) - q^B v (m 2 , q 2 ) . (A.4) 

Explicit calculation shows that = \T^ V (to all orders in (d — 4)). Defining B(m 2 ,q 2 ) by 



B(m 1 1 q 2 ) = B(m 2 , 0) + B(m 2 , q 2 ^ 

one finds for the relevant independent integrals 

A{m 2 ) 



(A.5) 



B(m 2 ,0) 



16vr 2 



(A.6) 



B(m 2 ,q 2 ^ 



I _IJl_4m£ l g ( ]±V}E^R 



~ifev 1 _ 4 ™7<? 2 % 2 - 4 ™ 2 



8tt- 



4m 2 



1 tan 



-i 



4m 2 



-1 



1 — 4m 2 jq 2 tanh 



-i 



■y/l-4m 2 /q 2 



(g 2 > 4m 2 ) 

(0 < q 2 < 4m 2 ) 

(g 2 < 0) 



(A.7) 
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and, with T^m 2 , q 2 ) = Ti(m 2 , q 2 )(q^q v - q 2 g VLV ) + T 2 (m 2 , q 2 )g^ u , 



Ti(m 2 ,q 2 ) 



4B 21 (m 2 ,q 2 ) + ^p- + 0(d-4) 
m 2 



(A.8) 



T 2 (m 2 ,g 2 ) =2A{m 2 ) (A.9) 

with B 2 i(m 2 ,q 2 ) as defined in the text. 

The insertion graphs of Fig. 2 involve the integrals C(m 2 , q 2 ), C^(m 2 , q 2 ) and C^ u (m 2 , q 2 ) 
defined by 

i r r r 1= f ^ k {1, k^k v } 

1 ' (27r) rf p 2 -m 2 ) 2 ][(A;-g) 2 -m 2 ] ' 1 } 

These integrals occur in the calculation only in the combination 

U^ u = AC^ - 2q fl C u - 2q u C fl + q^q u C 

= Ui(m 2 , q 2 )^^ - q 2 g ilv ) + U 2 (m 2 , q 2 )g^ , (A.ll) 

where one may show that 

U x {m\ Q 2 ) = - B(m ^ q2) + 0(d - 4) (A.12) 
U 2 (m 2 , q 2 ) = - -±- + G(d - 4) . (A.13) 

As explained in the text, the 0{d — 4) terms in the expressions for the various integrals 
do not enter the final result and hence are not displayed explicitly. 

Finally, in order to recast the one-loop result in terms of the physical kaon masses, one 
requires the value of the derivative of Ti(m 2 , q 2 ) with respect to m 2 , which is readily obtained 
from (A.8), (A.l) and the relation 

d s / 2 2\ B(m 2 , q 2 ) i ,»-.,.>> 

; B 21 (m 2 , q 2 ) = + — — — . (A.14) 



dm 2 Z1V ^ ' 2q 2 192vr 2 m 2 
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FIGURES 
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(c) 

FIG. 1. One-loop contributions to the correlator II?,* 




(b) 

FIG. 2. Non-contact "insertion" graphs of 0(q 6 ) 
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(b) 

FIG. 5. Contact, non- "insertion" graphs of 0(q 6 ) 



FIG. 6. Tree-level graphs of 0(q 6 ) 
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FIG. 7. Loop contributions to the correlator II 38 (g 2 ) in units of 10~ 5 . The solid line is the 
one- loop result, the dotted and dashed-dotted lines the imaginary and real parts of the full 
one- loop-plus- two-loop result, respectively. 
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